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Abstract. This paper reviews quantum oscillation studies of iron-pnictide 
superconductors and related materials. These measurements give unique information 
regarding the full three dimensional topology of the Fermi surfaces and the 
renormalisation of the quasi-particle masses. The review will cover measurements 
of the 122 arsenide end members, XFe2As2 (X = Ba, Sr, Ca) which have a spin density 
wave ground state, but will concentrate on the phosphide end members (LaPePO 
and XPe2P2) which have a paramagnetic ground state and a Permi surface topology 
which is similar to the higher Tc superconducting iron-pnictides. All three of the 
122 phosphides become superconducting when P is partially substituted by As, and 
for the BaFe2(Asi_a;Pa;)2 series de Haas- van Alphen oscillations are observable for 
0.42 < a; < 1 with Tc up to 25 K. The results show the changes in the Permi surface 
topology and the increase in the mass renormalisation as the correlations which induce 
superconductivity develop. 



1. Introduction 

A popular view of tlie iron-pnictide superconductors is that their unusual properties 
stem from the unique structure of their Fermi surfaces. Very soon after the discovery 
of superconductivity at temperatures greater than 40 K, conventional density functional 
theory (DFT) band-structure calculations using the local density approximation (LDA) 
(and its derivatives such as the generalised gradient approximation (GGA)) showed 
that the Fermi surface is mostly composed of small, almost two dimensional, tubes 
running along the c-axis [1]. The tubes at the centre of the Brillouin zone are hole-like 
whereas those at the corner of the zone are electron-like [see Figs. ([1]), (|5]) and ([7])]. 
For the undoped parent materials, for example LaFePO or LaFeAsO, the volume of the 
hole and electron Fermi surfaces are exactly the same and the material is said to be 
compensated as the number of electrons and holes are equal. 

It was noticed that as these Fermi surfaces are both roughly circular in cross section 
and as they are the same size, then a single translational vector Q = (tt/o, vr/a) maps one 
surface onto the other. This property, known as nesting, causes a peak in the Lindhard 
dielectric response function x{q) q = Q which can lead to charge or spin density 
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wave instabilities depending on the relative strength of the electron-phonon or electron- 
electron interactions. If the nesting is not so perfect, and assuming electron-electron 
interactions dominate (as they do for iron-pnictides) then strong spin fluctuations should 
be expected with wavevectors centred around Q. In many models [21 [Sj HI [5] it is these 
spin-fluctuations which provide the pairing interaction which binds the electrons into 
Cooper pairs and causes superconductivity. As the Fermi surface geometry is at the 
heart of this interaction its experimental determination is of high importance. 

It should be said before continuing that a completely itinerant electron model, as 
outlined above, is unlikely to be the whole picture. Although the many body electron- 
electron correlations are usually assumed to be relatively weak in the iron-pnictides, as 
compared to the cuprate high superconductors for example, they cannot be neglected. 
For example, according to the DFT calculations the Fermi surfaces of LaFePO and 
LaFeAsO are very similar, however experimentally LaFePO is a paramagnetic metal 
which becomes superconducting below 6 K [6] whereas LaFeAsO has a transition to spin 
density wave metallic antiferromagnetic state below 130 K and is not superconducting 
[7]. An explanation for this difference is to suppose that the magnetic interactions 
originate primarily from many body electron-electron correlations (for example the U 
and J parameters in Hubbard type models) but the structure of the spin-fluctuations 
in momentum space is governed by the Fermi surface geometry [8]. This is backed up 
by the fact that DFT models generally correctly predict the form of the magnetic order 
in the static SDW phase. The relative importance of the local interactions and the 
itinerant behaviour is still under debate. 

Given the importance of determining the properties of the Fermi surface it is 
not surprising that this has been the subject of many theoretical and experimental 
investigations. By far the most popular and ubiquitous experimental technique is 
angle resolved photoemission spectroscopy (ARPES). Here the sample is illuminated 
with photons, typically of energy 6 — 100 eV, and emitted photoelectrons are collected 
and analysed. As the in-plane momentum of the electrons is conserved when they are 
emitted, measurements of the energy and momentum of these electrons as a function of 
angle can be used to deduce the e{k) relationships for the various bands in the metal. 
The technique yields a direct 'image' of a 2 dimensional slice through the Fermi surface 
at a particular kz which is determined by the energy of the photons (with an uncertainty 
determined by the escape depth of the photoelectrons). One problem with the technique 
is that it is surface sensitive. The surface electronic structure may differ from the bulk 
for several reasons. For example, the surface may be charged (and hence effectively 
doped) if the cleavage plane is not neutral. The technique also has significantly lower 
resolution than quantum oscillations studies (particularly in k^)- 

1.1. Determining the Fermi surface via quantum oscillations 

Measurements of the Fermi surface properties based on magneto-quantum oscillation 
(QO) effects are an important complementary technique to ARPES. Although the 
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majority of the Fermi surface parameters which may be determined by the QO effects 
are also accessible by ARPES, the main advantages of QO are: i) it is very insensitive 
to surface states as it is dominated by the bulk response thus giving an unambiguous 
probe of the bulk Fermi surface, ii) it has very high fc-space resolution in all three 
crystallographic directions - typically Fermi surface cross-sections can be determined to 
better than 10~^ of the area of the Brillouin zone, iii) it has very high energy resolution, 
so that it gives the effective mass very close to the Fermi level {AE ~ =0.6 meV 
for B—IOT). A major limiting factor is that only samples with sufficiently long mean 
free path will give a measurable signal. High magnetic field and low temperatures are 
usually essential. 

The technique is based on the fact that in a magnetic field the only allowed values 



the single electron flux quantum {h/e) and n is an integer known as the Landau index. 
These two dimensional electron tubes run parallel to the field B and for almost free 
electrons will have a circular cross-section. As the field is varied the k space radius of 
a particular Landau level n will increase, and as each level passes through the Fermi 
surface the states will depopulate giving rise to a saw-tooth like variation of the free 
energy, periodic in inverse field 1/B. The frequency (in inverse field) of the oscillations 
is related to the fc-space cross sectional area of the Fermi surface by the Onsager 
relation, F = hAk/2Tie. For a three dimensional metal will vary as a function 
of kz and hence each kz slice (assuming B\\z) will give rise to oscillations of varying 
frequency and phase. The total signal will be the sum of the contributions from each 
slice. Usually the sum is dominated by the extremal cross- sections of the Fermi surface 
where dA^/dkz — and hence a typical Fermi surface sheet will give rise to a small 
number of oscillatory signals. 

The oscillations in the free energy may be determined by measuring many different 
physical properties such as the resistivity (Shubnikov-de Haas effect), the magnetisation 
or magnetic torque (de Haas- van Alphen effect), specific heat etc. Finite temperature, 
impurity scattering and other effects smear out the discontinuity as the Landau levels 
pass through The response may be decomposed into a Fourier scries, and the signal is 
usually dominated by the fundamental harmonic at frequency F. For a three dimensional 
metal of volume V the full expression for the first harmonic of the oscillatory torque T 
is (in SI units) 
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This is commonly known as the Lifshitz-Kosevich (LK) formula [9]. The expression 
includes the most common attenuation factors, Rt, Rd and Rs which are due to finite 
temperature, impurity scattering and spin-splitting respectively. It should be noted that 
there are other attenuation factors, due to effects such as mosaic crystal structure or 
superconductivity which are important in some circumstances. In the expression the 
quasiparticle mass occurs several times in different forms, mg is the free electron mass 
and rrLh is the band mass as calculated from mean-field DFT band structure calculations. 
The mass m* which occurs in the expression for the temperature dependent damping 
factor Rt is the quasiparticle effective mass renormalized by electron-phonon and 
electron-electron interactions, rris is the 'spin-mass' which is renormalized by electron- 
electron interactions (which enhance the Pauli susceptibility) only. The curvature factor, 
(|^)~^/^ accounts for the number of states which contribute to the extremal orbit, here 
is the component of k parallel to the field. If changes little with fcy, for example 
because the Fermi surface is locally almost two dimensional, then these orbits will have 
a large intrinsic amplitude. 

The Dingle factor, Rd, is determined by the impurity scattering rate and is 
not normally renormalized by many body factors. For free electrons we can make the 
substitutions, mbVp = hkp, and i = tvf, where vp is the Fermi velocity and i is the 
mean free path, so that 

„ / ■nhkp\ , . 

For the general case, £ and kp should be replaced by values averaged over the orbit. A 
simple approximation is to take 7r/c|, = A = 2T:eF/h, so that Rd = exp(— 1140V^ / d-B), 
with I in Angstroms and F, B in Tesla. 

This paper will review QO experimental studies of the Fermi surface of the iron- 
pnictide superconductors. Because of the need for high purity samples these studies have 
mostly been conducted on the stoichiometric end members of the various iron-pnictide 
series. Charge doping of these parent materials, which is a common route to induce 
superconductivity, generally induces significant disorder which shortens the mean-free- 
path to such an extent that QO are not observable even in the highest available fields 
(-Bmax — TOT). On the other hand, isovalent substitution of As by P in the series 
BaFe2(Asi_a;Pa;)2 does not induce so much disorder and QO are observed for dopings 
producing up to ~80% of the maximum [10]. 

2. Parent iron-arsenides in SDW state 

In their undoped state, at ambient pressure, the '1111' (e.g., LaFeAsO) and '122' (e.g., 
BaFe2As2) iron-arsenides have a spin-density wave ground state. Below a temperature 
of 140 K BaFe2As2 undergoes a structural and magnetic transition (at slightly different 
temperatures) to a state with stripe ordered magnetic moments with a unit cell 
which is roughly a/2 larger (along the edge) and rotated by 45° compared to the 
high temperature cell [llj. The structure is actually orthorhombic (Fmmm), with 
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Figure 1. Simplified scliematic picture of Fermi surface reconstruction due to band 
folding at the spin-density wave transition in the iron-pnictides. a) Shows circular 
Fermi surface sheets, hole like at the zone centre (F) and electron like at the zone 
corner (M), in the paramagnetic Brillouin zone (solid line). The change in structure 
at the SDW transition will cause the Brillouin zone to change (dashed line), b) shows 
the band-structure energy - momentum curves corresponding to the Fermi surface in 
a). As the Brillouin zone changes the band at M is translated (folded) back to the F 
point in the reduced zone scheme, c) The folded bands hybridize and a gap forms at 
the Fermi level ep, and the Fermi surface disappears. 



a=5.6146A, 6=5.5742 A,and c=12.9453A. The fact that a and b only differ by 0.7% 
means that the crystals are highly susceptible to twinning, so that crystals usually have 
multiple domains in which the a and b directions are switched. Twin free samples have 
been produced by cooling through the SDW transition with an applied uniaxial stress 
or a large magnetic field [121 1131 [15 • 

As a first approximation the effect of this change in structure on the Fermi surface 
can be understood by a simple band-folding argument. As illustrated in Fig. [H the 
change in structure causes the electron band which produced a Fermi surface at the 
corner of the original cell to fold to the centre of the zone. The electron-like bands and 
the hole-like bands now cross and in the simple case where the shape and size of hole 
and electron Fermi surfaces are identical, the bands would cross exactly at the Fermi 
level. The bands will then hybridize and a gap will appear at the Fermi level, hence the 
Fermi surface would completely disappear and the material would become insulating. 
In reality, the original Fermi surface sheets do not have exactly the same in-plane shape 
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Figure 2. Quantum oscillations in BaFe2As2. a) Raw oscillator frequency shift which 
is proportional to the resistance versus field for B\\c and B J- c (upper and lower curves 
respectively), b) Fourier transform of the i3||c data, c) Variation of observed SdH peaks 
with field angle (solid lines are fits to a simple elliptical Fermi surface model), d) DFT 
fixed moment calculation of the Fermi surface using the negative U procedure to fix 
the ordered moment. The predicted extremal orbits are labelled 1-4. Figure adapted 
from Analytis et aL[16j. 



and there is some corrugation along k^, hence the Fermi surface will not be completely 
wiped out and small three dimensional Fermi pockets will remain. In addition, the 
formation of the SDW state causes the bands to shift in energy (some bands more than 
others) so the exact band-structure is not simply related to the folded paramagnetic 
structure. Fixed moment DFT calculations seem to produce approximately the correct 
structure - although not the exact Fermi surface topology [15]. 

Quantum oscillation experiments have been reported on several members of the 
'122' family XFcgAss with X=Sr, Ba and Ca [HI [IH [H]. Examples of the raw 
(background subtracted) experimental results for the oscillations in the resistivity of 
BaFe2As2 (measured inductively using a tunnel diode oscillator) are shown in Fig. |2l The 
signals are relatively weak and can only be observed in very high magnetic field (typically 
greater than SOT). In all these measurements only very low oscillation frequencies were 
observed. For example for BaFe2As2, the values with i?||c were reported to be 80 T, 190 T 
and 440 T, which correspond to 0.3, 0.7 and 1.7% of the total area of the paramagnetic 
Brillouin zone respectively [19]. The weak variation of the a frequency with magnetic field 
angle shows that this comes from a small quasi spherical pocket as anticipated. A DFT 
calculation of the Fermi surface of BaFe2As2 based on the observed magnetically ordered 
unit cell is shown in Fig. These calculations use the LDA+U procedure, which is 
common for strongly correlated metals, with a negative orbital potential U to reduce the 
ordered moment. This negative U has no physical significance and should be regarded as 
an empirical fix to the LDA calculations. With moderate shifts of the calculated bands 
energies (< 60meV), the predicted orbits 1-3 can be brought into agreement with the 
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Figure 3. Quantum oscillations in LaFePO. a) Raw torque versus field at T=0.35K 
with 4*'' order polynomial subtracted, b) Fast Fourier transform of the same data, 
c) FFT spectra of different samples (labelled A, L D) at several angles 9 showing the 
various extremal orbit peaks. Figure adapted from data in Ref. |21j . 



experimental results in Ref. [16]. In this study the Fermi surface determination was 
incomplete as all the predicted orbits were not observed experimentally. 

Recently a complete Fermi surface determination of de-twinned BaFe2As2 was 
reported |20j. In these detwinned samples the oscillations are large in amplitude even 
in low fields less than 10 T. All extremal orbits of the calculated Fermi surface were 
observed, and the calculations could be brought into almost exact agreement with 
experiment with appropriate band energy shifts (< 65meV). The determined Fermi- 
surface is quite similar that shown in Fig. |2] but has some differences in detail for some 
of the pockets. 

3. LaFePO 

The 1111 compound LaFePO was the first discovered iron-pnictide superconductor |6] 
and was also the first paramagnetic iron-pnictide in which QO were observed [21]. Unlike 
the isomorphic arsenide LaFeAsO, LaFePO is superconducting in its stoichiometric state 
(some reports have suggested that a few percent of oxygen vacancies are needed to induce 
superconductivity [22l |23], however, it is not clear if this is necessary [21]). Unlike 
LaFeAsO, LaFePO is not magnetically ordered and magnetism can not be induced by 
small amounts of doping away from stoichiometry. The superconducting transition 
temperature is quite low Tc ~ 6K and the upper critical fields {Hc2\\c = 0.68 T, 
i?c2 -L c = 7.2 T) are easily exceeded with conventional superconducting magnets [23]. 
It is therefore an ideal candidate for study by quantum oscillations effects. 

Data were first reported by Coldea et al [21], closely followed by Sugawara et 
a/. [25] As the best available single crystals of LaFePO are quite small (typically 
planar dimensions of a few hundred microns and thickness a few tens of microns) the 




Figure 4. Observed dHvA oscillation frequencies in LaFePO (solid symbols) versus 
field angle as the field is rotated from B\\c {9 = 0°) to B _L c. The values derived 
from the DFT calculation, after shifting the band energies (as described in the text) 
are shown by the solid lines. The right panel shows the extremal orbits on the electron 
and hole Fermi surface sheets. Figure adapted from data in Carrington et al [241 . 



microcantilever torque method is well suited. Examples of raw torque versus magnetic 
field data for this compound are shown in Fig. [31 The data are best appreciated by 
calculating the fast Fourier transform (with respect to inverse field). A total of 7 
fundamental frequencies are observed plus harmonics of the strongest amplitude signals 
[a and /3 orbits). Measurements were repeated as a function of field angle as the field 
was rotated from parallel to the c-axis {6 = 0°) towards the a6-plane. As the field is 
rotated the extremal electron orbits transverse different sections of the Fermi surface, 
and hence by analyzing the field angle dependence of the dHvA frequencies, F{6), a 
detailed model of the Fermi surface can be constructed. 

For a perfectly two dimensional Fermi surface, F is inversely proportional to cos 6' 
and hence by plotting F cos 6 versus 6 small deviations from two dimensionality can be 
easily seen. The simplest distortion of a two dimensional Fermi surface in a tetragonal 
system is a cosine dispersion, where the e{k) = fi^{k1 + ky)/2m + t^cos(cA;2). In this 
case, there are two extremal frequencies, the minimum -Fmin and the maximum -Fmax, 
which vary with angle according to the Yamaji formula |26] 



mm, max 



(6*) cos 6* = Fav±AFJo{ckFtan( 



In this expression, the first term is the mean frequency {Fav = -Fmin + -^max)/2 at 6* = 0, 
the prefactor of the second term is AF = {F^^^ — F^jjj)/2, Jq is the Bessel function 
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Figure 5. Calculated electronic structure of LaFePO without spin orbit interaction 
or adjustment to fit experimental data, a) Band structure energy-momentum along 
the zone diagonal. The percentage of Fe-d character is indicated by the colour scale, 
with the band with predominant Fe (1^2 character overprinted in black, b) Brillouin 
zone with symmetry labels, c) Three dimensional Fermi surface, d) c axis (110) slice 
through the Fermi surface with predominant band character indicated by the colour 
scale. Figure adapted from Ref. [24] . 



and kp = 2eFav/h. At certain angles Jo = and F^i^ = Fmax and here all parts of the 
Fermi surface interfere constructively to give a large increase in the amplitude of the 
oscillations. Such a point is observed at 6* = 42(2)° for the (3 frequencies of LaFePO, 
which together with the overall behaviour of F^i^ and Fmax are very strong evidence that 
they originate from a single c-axis warped section of Fermi surface (see Fig. H]). The two 
a frequencies are much closer indicating that they come from an almost two dimensional 
section of Fermi surface with little warping, and no clear Yamaji amplitude peak is seen. 
For the lowest frequency 6, Fcos6 shows almost no change with 6, consistent with it 
coming from an almost 2D Fermi surface section. The remaining two frequencies e and 
7 are much weaker; they were only observed above 5=30 T and again Fcos9 shows 
minimal change with 6 within the noise. 
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Figure 6. In plane cross sections of the calculated Fermi surface of LaFePO including 
the spin-orbit interaction and with bands shifted to best match the dHvA data. Figure 
adapted from Ref. [24]. 

Further analysis of the data is greatly aided by band-structure calculations. 
The experimental dHvA frequency versus field angle data can be compared to these 
calculations and adjustments made to best fit the data. Often this can be accomplished 
with a small number of parameters, usually a single energy shift for each band or group of 
bands. The alternative is to fit the data to a series of appropriate geometrical harmonics. 
This latter approach works well when the Fermi surface is relatively simple and only 
has small deviations from a sphere or a two dimensional tube and only a small number 
of harmonics are required (for example, SrRu204 [27j or Tl2Ba2Cu06+5 [28| 129]). 

For LaFePO the band-structure calculations were done [211 El] using the augmented 
plane wave plus local orbitals methods as implemented in the WIEN2K code [30]. For 
these calculations the experimental lattice parameters and internal coordinates were 
used. A plot of the energy- momentum band structure along the T — M direction (Fig. [5]) 
shows two hole bands crossing the Fermi level close to F and two electron bands crossing 
close to M. So apart from the doubling of bands the detailed band-structure looks quite 
similar to the simplified picture in Fig. [1] There are however two complications. First 
the outer electron surface has significant c-axis and in-plane warping and second, there is 
a fifth band, mostly coming from the Fe dz^ orbital which cuts across the two dimensional 
hole bands and forms a closed three dimensional pocket close to Z. 

The assignment of the a and /3 orbits to the inner and outer electron sheets can be 
made rather easily. The calculation fairly accurately predicts the size and warping of 
these sheets. The relatively strong c-axis warping of the /3 electron sheet and the much 
more two dimensional nature of the a electron sheet is well reproduced. To get exact 
agreement with the data it is necessary to shift the calculated electron bands up in energy 
(-1-85 meV and -|-30meV for the a and /3 orbits respectively). The physical meaning of 
these band shifts will be discussed later. With this adjustment the agreement is very 
good. The calculated c-axis warping of the /3 sheet is slightly too large: experimentally 
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AF/ Fav=0.23 whereas the calculation gives AF/Fav=0.33, but for the a sheet it is 
much closer. Without the inclusion of spin-orbit coupling the two electron bands are 
degenerate along the zone edges X — M and R — A. This degeneracy is lifted by the 
spin-orbit interaction and a gap of ~ 50meV appears in the calculations. As the two 
electron bands need to be shifted by different amounts to get agreement with experiment 
this implies that this gap is actually about double this (i.e., around 100 me V). This gap 
can be seen on the cross sections through the fitted Fermi surfaces shown in Fig. [61 

Identification of the hole orbits is difficult. The lowest observed frequency S is 
likely the minimum of the inner hole sheet which is centred at F. As predicted by the 
calculation there is minimal change in F cos 6 for this orbit as the sheet is locally close 
to two dimensional at this point. It is brought into perfect alignment by shifting all the 
hole bands down by 53meV. In that case the maximum of the inner hole band almost 
perfectly coincides with the inner electron band a. The assignment of the other two 
frequencies 7 and e is more problematic. With the above mentioned shift of the three 
hole band energies it can be seen (Fig. H]) that 7 and e approximately correspond to the 
minimum and maximum of the largest hole sheet. The variation of e with 6 however 
is much less than expected. As these orbits are only observed at high magnetic field 
(>30T) it is possible that either or both of them result from magnetic breakdown or 
torque interaction effects arising from the electron orbits [21]. In the calculations, the 
energy of the Fe dz^ band, which gives rise to the warping of the hole sheets and the 
small spherical pockets at Z, varies very strongly with the height of the P atom and so 
it is possible that it does not actually cross the Fermi level. In this case the above rigid 
band shifting will not work accurately because the d^^ band has crossed and hybridized 
with the two dimensional dxz\\xy bands. Instead, we might expect the hole bands are 
almost perfectly two dimensional, thus giving rise to only two dHvA frequencies with 
no variation in FcosO as a function of 6. 6 is almost certainly from the inner hole band 
but which of 7 and e is from the largest hole band is uncertain. 

The effective mass of the quasiparticles on each of the orbits was deduced by 
measuring the temperature dependence of the dHvA oscillation amplitudes and fitting 
this to the expression for (Eq. (1)). The masses for the electron orbits range from 
1.7 to 1.9 me- This should be compared to the masses calculated from the band- 
structure (for the shifted bands) which range from 0.7 to 0.9 mg, hence the many-body 
mass enhancement factors m*/mh — 1 = A range from 1.0(2) to 1.7(2) for the 4 orbits 
(see Table [1]). The correlation is much stronger than expected from electron-phonon 
scattering Agp — 0.21 |3T], and likely results from electron-electron interactions. ARPES 
measurements on LaFePO concluded that the bands were all narrowed by a factor ~ 2.2, 
hence this implies that the electron correlations act over a wide energy range rather than 
being localised close to the Fermi level. This goes against interactions with low energy 
spin fluctuations being the major cause of the mass enhancement. The masses of the hole 
orbits are quite similar, however the errors are larger on account of the smallness of the 
signals. It is unclear at present whether the masses of the observed orbits are sufficient 
to account for the measured electron specific heat 7e. Four two dimensional Fermi 
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Figure 7. Fermi surfaces of the 122 iron phosphides, BaFe2P2, SrFe2P2, and 
CaFe2P2, determined by DFT calculations with bands shifted to best match the dHvA 
measurements [36 l [37 l [38 ] , 

surfaces each with m* /nte = 2 would give 7e = 6mJ mol~^ K^^ whereas the reported 
value on polycrystalline samples was 7e = 12.7 mJ mol^^ K^^ [22j. This discrepancy 
could indicate that a heavy hole sheet {m* /rrie — 11) has not yet been detected but the 
heat capacity measurement also needs to be checked on single crystal samples. 



4. 122 phosphides: BaFesPa, SrFeaPa, CaFesPa 

Studies of the Fermi surfaces of the 122 phosphides are important for a number of 
reasons. Unlike the corresponding arsenides, the Fermi surfaces of the phosphides are 
very similar to those of the superconducting phases of the Ba and Sr 122 materials. In all 
three (Ba, Sr, Ca) phosphides superconductivity can be induced by partial substitution 
of P with isovalent As [321 [33] and hence the pure phosphide can be viewed as the 
end member of these series. The maximum produced by P/As substitution is very 
comparable to that obtained by the non-isovalent dopants K fTT] / Co fM] or pressure 

The 122 phosphides are all paramagnetic, non-superconducting metals. Calcula- 
tions of the Fermi surface of the Ba and Sr phosphides closely resemble those of the 
corresponding arsenides in their paramagnetic state (i.e., above Tsdw)- The Fermi sur- 
face of CaFeaPa on the other hand, closely resembles that calculated for CaFeaAsa in its 
high pressure state. One effect of replacing As by P in the 122 materials is to decrease 
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Figure 8. Experimental dHvA frequencies versus magnetic field angle for rotations 
from B\\c [001] to B\\a [100] for both SrFe2P2 [36 and BaFe2P2 EE] (symbols). The 
unlabelled points in grey in the BaFe2P2 panel are harmonics of the a and /? orbits. 
The orbit labels for SrFe2P2 have been changed from that in Ref. j36] to be consistent 
with those used for BaFe2P2- The solid lines are the frequencies calculated from the 
DFT calculations with shifted bands. The far right panel depicts the locations of the 
various orbits on the Fermi surface of BaFe2P2- The labels also apply to SrFe2P2 which 
has a very similar Fermi surface (see Fig. [7]). Figure adapted from Refs. jSGj and [3S] 

the distance between the Fe planes and hence increase the pnictogen-pnictogen bonding. 
In the case of CaFe2As2 the apphcation of pressure also decreases the interplane distance 
and for P > 0.35 GPa there is sudden, very large (~10%) decrease in the c-axis length. 
This high pressure phase is known as the collapsed tetragonal state. CaFe2P2 has a 
similar interplane distance as the collapsed tetragonal state of CaFe2As2 and hence a 
similar Fermi surface. The calculated Fermi surfaces for all three compounds are shown 
in Fig. [3 

dHvA studies have been completed for all three of these 122 phosphides [361 EZl 
[TOl l38] . The oscillations can be observed down to low field (~ 5T) and clear dHvA 
peaks in the Fourier spectra are observed for orbits on all four sheets of Fermi surface. 
Experimental dHvA frequencies versus magnetic field angle for both SrFe2P2 [36] and 
BaFe2P2 [M] are shown in Fig. |H1 For SrFe2P2 by far the largest signal comes from the 
a (inner electron) orbits, followed by the /3 (outer electron) orbits. The large a orbit 
signal originates from the fact that the in-plane area of this orbit changes very little 
with kz and hence the curvature factor (Eq. |2]) is large. This is demonstrated by the flat 
behaviour of Fa cos{9) in Fig. [HI It should be noted however, that in this case, this does 
not mean that this sheet of Fermi surface is highly two dimensional. Because of the screw 
symmetry of the body centred tetragonal structure of the 122 materials, the in-plane 
elliptical shape of the Fermi surface at the corner of the zone rotates by 7r/2 between 
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Table 1. Mass enhancement factors (A — m* /mi, — 1) in iron-phosphides. Here m* is 
the measured thermal mass and mi, is the DFT band structure mass at the same angle 
for with the band shifted to best match the observed dHvA frequency. For Ba-122, 
Sr-122 and LaFePO the orbit labels are as in Figs. ^ and (g]). For Ca-122 the hole 
orbit labelled here as 62 is the /3 orbit in Fig. (fTO|) . 
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the centre and top of the zone (see Fig. [9]). This type of distortion has been described 
as resembling a snake that has swallowed a chain [2^. Exactly half way up the zone, 
this odd symmetry cosc/c^sin20 term (/c2,i in the notation of a cylindrical harmonic 
expansion [271 [29] ) is zero, and the shape is closest to circular (even symmetry cosine 
distortion terms remain) (see Fig. [9]). The anisotropy of the individual Fermi surface 
sheets, parameterized by the anisotropy of the plasma frequencies (which in the isotropic 
scattering rate limit is equal to the square-root of the resistance anisotropy) are shown 
in Table H 

For SrFe2P2 the calculated (a//3) electron Fermi surfaces are brought into almost 
perfect alignment with the data by shifting them up by 49 meV and 59 meV respectively. 
The outer hole sheet in this material has very strong warping close to the top of the 
zone which as in LaFePO is caused by the band acquiring strong Fe character. At an 
angle of ^ ~ 50° the large warped sections of Fermi surface are almost parallel and so 
there the curvature factor strongly increases. Hence, it is close to this angle where this 
orbit is most easily observed. Interestingly, the calculation correctly predicts the size of 
this part of this Fermi surface without any band energy shifting necessary. As will be 
discussed below this likely is because this Fermi surface section is does not nest with 
any of the electron sheets. The inner hole sheet is somewhat smaller than the electron 
sheets and is shifted down in energy by 110 meV compared to the DFT calculation. This 
causes the sheet to pinch off (see Fig. [7]), so 5i in Fig. [8] ceases to exist. 

An important check of the band shifting fitting of the DFT calculations to the 
dHvA data is to calculate the volume of the Fermi surface sheets. These materials are 
all expected to be compensated metals, so they should have equal numbers of electrons 
and holes. For SrFe2P2 with the shifts, the 6, 7, a and /3 bands contain -0.026, -0.284, 
0.108 and 0.197 hole/electrons respectively, and so the charge balance remains almost 
exact. The imbalance of 0.005 extra holes per unit cell is equivalent to < 2% of the 



Quantum oscillation studies of the Fermi surface of iron-pnictide superconductors 15 




010 



Figure 9. Cross-sections of the determined Fermi surface of BaFe2P2 in the ab plane 
at three different values (quoted in units of c*). The strongest band character at 
each fc-point is indicated. Note that the zone corner is labelled X in the body centred 
unit cell of the 122 compounds. This is equivalent to M in the simple tetragonal unit 
cell of LaFePO. Figure adapted from Ref. [38j . 

Table 2. Calculated in-plane (w* ,j) and c-axis {uj* J plasma frequencies in eV 
units using the bands shifted to best fit the dHvA data. The values have been 
renormalized by factor (m*/TOf,)^^/^ where m* /mt, is the average mass enhancement 
for the particular sheet (see Table [1]). The anisotropy F = uj*^/u!*^, which in the 
isotropic scattering approximation is the square root of the conductivity anisotropy, is 
also given for each Fermi surface sheet. For LaFePO the possible contribution of the 
small spherical sheet has been ignored. The anisotropy of the remaining two hole 
sheets is a lower bound as all orbits have not clearly yet been identified. 
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total volume of the hole sheets. 

The data for BaFe2P2 are in many ways quite similar to those for SrFe2P2- The 
general shape of all four sheets of Fermi surface are very similar and also the absolute 
sizes are quite comparable. The most significant difference is the size of the inner hole 
sheet, particularly orbit 62- The cross section of this sheet is almost identical to the a 
electron sheet. In fact, the frequency 62 is right between the two extremities of that 
from the a sheet. This suggest the possibility that there could be almost perfect nesting 
between these electron and hole sheets. As mentioned above, however the electron 
sheets have quite elliptical cross-sections at the centre and top of the zone so although 
the areas match the shape does not. However, half way between these two extremities 
the electron sheet becomes almost circular and therefore matches very well the hole 
sheet (see Fig. [9]). In addition the band-character of these sheets is also very similar. It 
should be recalled that despite such good nesting BaFe2P2 is neither superconducting 
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CaFe2P2 




G (degree) 

Figure 10. Experimental dHvA frequencies versus magnetic field angle for rotations 
from B\\c (001) to B\\a (100) for CaFe2P2 ED- The electron and hole Fermi surface 
are shown to the right of the figure. Figure adapted from Ref. [37] . 

nor magnetic. Hence this is a good example showing that nesting is a necessary but not 
a sufficient condition for the formation of those ground states. Many body correlation 
effects which become stronger for the arsenide materials clearly also play a key role. 

The band shifts needed to get the DFT calculation to agree with the data for 
BaFe2P2 (+68/+58meV for the inner/outer electron sheets and -113 meV for the inner 
hole sheet) are almost identical to those for SrFe2P2- In particular, again it is found that 
the highly warped d^^ section of the outer hole sheet agrees well with the calculation 
without any shift needed. 

Next we turn to CaFe2P2- Here the dHvA spectra are significantly simpler than the 
other 122 compounds with far fewer observed frequencies (see Fig. [T0|) . As remarked 
above (Fig. [71), the Fermi surface of this material is markedly different from the other 
122s. The hole Fermi surface has changed from two quasi two dimensional cylinders to 
a single flat pancake which is depressed close to the T ~ Z line. The electron sheets 
have also become much more warped but still roughly keep their cylindrical shape. For 
this material, it is found that the DFT calculations agree almost perfectly with the 
dHvA data without any adjustment of the band energies (to with ± lOmeV). This is an 
important counter-example to the results found for LaFePO, BaFe2P2 and SrFe2P2- It 
would therefore appear to be clear that the band shifts are not an inherent inaccuracy 
in the mean-field band-structure but are rather caused by non-mean-field correlation 
effects which act between almost nested Fermi surfaces. A obvious candidate for this 
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are spin-fluctuations. 

The mass enhancements found from the dHvA measurements for all three 122s 
plus LaFePO are shown in Table [H CaFe2P2 has the smallest enhancements, which are 
about two times larger than those expected from conventional electron-phonon coupling 
{Xep — 0.23) [39]. For SrFe2P2 the enhancements are quite sheet dependent, being 
only slightly larger than for the Ca compound for most sheets except for the most two 
dimensional a electron sheets where they are a factor of two larger. BaFe2P2 again has 
some sheet dependence and the average A = 0.74 - about 50% larger than for CaFe2P2- 
A useful check of consistency is to apply the measured mass enhancements for each 
sheet to the calculated density of states for each sheet and therefore get an estimate of 
the electronic specific heat. For SrFe2P2 this gives je = 10.4(2) mJ mol~^K~^, which 
agrees well with the experimental value of 7^; = 11.6(2) mJ mol~^K~^ [40] . 

4.1. BaFe2(Asi.,Pj2 

Although studies of the non-superconducting end members described above can give 
valuable insight into the bulk electronic structure of the superconductors, it is clearly 
preferable to measure the actual superconducting materials themselves. Unfortunately, 
in the main, the condition that the samples must be free of disorder means that it is not 
possible to observe quantum oscillations in heavily doped superconducting samples. 
This is compounded by the fact that often the highest samples also have very 
high Hc2 values and so the normal state cannot be accessed. LaFePO and KFe2As2 
[IT] are rare examples where superconductivity occurs in the (close to) stoichiometric 
composition and dHvA oscillations are observable. However, it is unclear whether these 
low Tc materials are representative of the higher pnictides. Fortunately, there is one 
exception to this. In the isovalently substituted BaFe2(Asi_2:Pa;)2 series the residual 
resistance, and hence scattering, remains low even as the maximum Tc is approached 
(T^'""^=30K for X = 0.33) [12]. In this material, and other pnictides, there is clear 
evidence that the transport properties develop 'non-Fermi liquid' properties as the 
maximum is approached. The resistivity follows a power law as a function of 
temperature, p (x where the exponent n evolves from its conventional Fermi liquid 
value of 2 for the overdoped phosphide to close to 1 at optimal doping. Such behaviour is 
strongly reminiscent of the high cuprate superconductors [l3] and heavy Fermions in 
the vicinity of a quantum critical point [B]. In BaFe2(Asi_i:P2,.)2 the highest Tc occurs 
approximately where the SDW order temperature extrapolates to zero and hence is 
likely strongly influenced by the presence of a quantum critical point [15]. 

Using pulsed magnetic fields up to 60 T dHvA oscillations were observed in the 
BaFe2(Asi_a;Pa;)2 series for x down to 0.42, where has reached 25 K which is ~ 80% 
of the maximum [TO]. The quite remarkable fact that oscillations could still be seen 
in such a heavily substituted composition implies that As/P cross substitution induces 
very little disorder on the conducting Fe plane. The observed frequencies were identified 
as coming from the electron sheets, and for the lower values of x (higher T^) only a single 
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Figure 11. Experimental dHvA frequencies for the a and /? electron orbits of 
BaFe2(Asi_j:Pj:)2 versus the phosphorous content x. Where more than one a/ (3 
frequency was observed the average frequency /mass is plotted here. The x dependence 
of the effective mass m* of the /3 orbits and Tc is also shown. Figure adapted from 
Ref. [ini [38] . The more accurate value of m* from Ref. [38] for a: = 1 is plotted in 
preference to that in Ref. [TO] . 



frequency from each electron band was discernable. In the higher Tc samples oscillations 
were only observable at the highest fields and hence the Fourier peaks are too broad to 
see the individual peaks. 

Two important results were found. First it was observed that both the a and 
(3 frequencies decrease markedly with decreasing x (see Fig. [H]). No such decrease is 
expected from the DFT band structure calculations. Second, the effective mass increases 
strongly as the maximum T^. is approached. It is likely that both these changes are 
caused by the increase in correlations (spin fluctuations) as the SDW ordered state is 
approached. It is not clear from the available data whether the increase in effective mass 
is a divergence associated with a possible quantum critical point or rather just a general 
increase with decreasing x. Unfortunately, no data is yet available for the values of x 
below the critical value. 

Recently Analytis et al. |16] reported measurements for the x = 0.63 with Tc=7K 
where orbits on the smaller hole sheet were observed. The frequency of these orbits 
were found to correspond to almost the same cross-sectional the smaller electron 

pocket, very similar to that found for the x = 1 composition described above. Hence it 



Quantum oscillation studies of the Fermi surface of iron-pnictide superconductors 19 

would appear that the relative size of the electron and hole sheets remain remarkably 
constant as a function of x. As the electron sheets shrink with decreasing x the hole 
bands shrink by a similar amount without any major change in topology. This is 
supported by recent ARPES results for the x = 0.38 composition [47J which show 
that in this close to optimal doped composition the shape of the Fermi surfaces closely 
resemble the end member BaFe2P2- 

A semi quantitative model of the electron / hole pocket shrinking of the quasi- 
nested electron-hole Fermi surfaces of iron-pnictides has been suggested by Ortenzi et 
al. IH]. Here scattering caused by coupling of the itinerant electrons to a bosonic spin 
fluctuation mode (assumed to be Einstein-like) causes both electron and hole pockets to 
shrink and also the effective mass (and T^) to increase. This was applied to a simplified 
model of the LaFePO Fermi surface and reasonable quantitative agreement was found 
between the magnitude of the energy shifts (which in this model determines the strength 
of the scattering potential) and the mass enhancement and even Tc. It should be noted 
however that taken at face value this model would predict low-energy correlation-driven 
renormalisation of the band energies (i.e., kinks in e{k) close to the Fermi level), whereas 
a comparison of the dHvA data to ARPES suggests the renormalisation occurs over a 
much wider energy range. It would however be intriguing to see if an extension of this 
model could explain the correlations between band energy shifts, effective masses and 
Tc in the BaFe2(Asi_2:P2;)2 series. This could prove to be a decisive test of the spin 
fluctuation pairing model in the iron-pnictides. 

5. Summary and Conclusions 

Quantum oscillation measurements have proven to be a very valuable probe of the 
bulk electronic structure of iron-based superconductors. Combined with first principle 
electron structure calculations a very detailed model of the Fermi surface and the quasi- 
particle dynamics can be obtained. This review has shown that density functional 
theory calculations provide a good starting point. Using fc-independent shifts to the 
band energies the band structure calculations can be fitted to the data using very few 
parameters. These shifts appear to originate from strong electron correlations and 
are present even for the non-superconducting, non-magnetic phosphides. Importantly, 
it is found that only those parts of the Fermi surface where there is electron-hole 
quasi-nesting are shifted in energy. For CaFe2P2 no shifting is necessary. In the 
BaFe2(Asi_2;Pa;)2 series the shifting increases as the SDW phase transition is approached 
and is accompanied by a strong increase in the effective mass. It is hoped that in the 
near future, more sophisticated first principle calculations which go beyond the mean- 
field level, such as dynamical mean field theory, may explain these results and lead to 
new insights into the physic of iron-based superconductors. 
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